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The transverse nuclear magnetic resonance (NMR) decays ol 129 Xe in polycrystalline xenon were 
recently shown to have a universal property: in the long-time regime these decays all converge 
to the same sinusoidally modulated exponential function irrespective of the initial transverse spin 
configuration prepared by a sequence of one or more radio frequency pulses. The present work 
constitutes a more comprehensive survey of this phenomenon. It examines transverse decays for 
several different isotopic concentrations of 129 Xe, employs additional pulse sequences, and performs 
similar measurements in a different material: 19 F in single-crystal and polycrystalline CaF2. We 
additionally verified the polycrystalline nature of our frozen xenon samples by X-ray diffraction 
measurements. With the possible exception of polycrystalline CaF2 where the observation of the 
long-time behavior is limited by the experimental resolution, all these systems display the long-time 
universal behavior characterized by particular values of the exponential decay coefficient and beat 
frequency that were unique for each lattice. This behavior has been theoretically predicted based 
on the notion of microscopic chaos. 



I. INTRODUCTION 



Nuclear magnetic resonance (NMR) has yielded 
valuable information to researchers for over five 
decades. Motivations for manipulating the magnetic 
nuclear states of atoms range widely and cover a vast 
array of disciplines and applications. In solid-state 
NMR, spin-spin relaxation is one of the important 
sources of microscopic information. The disquieting 
fact is, however, that despite the ubiquitous applica- 
tions of the methods and results of solid-state NMR, 
the accuracy of first-principle calculations of many 
spin-spin relaxation related quantities continues to 
be limited. ^ While many theories have been pro- 
posed to describe NMR spin-spin relaxation under 
various manipulations, none has succeeded in pro- 
viding a controllable method for predicting measure- 
ments from a knowledge of the microscopic Hamil- 
tonian of the many-spin system. The properties of 
the spin-spin relaxation at times long compared to 
the characteristic spin-spin correlation time Ti (the 
"long-time" regime) are of particular interest in this 
respect because they are the most difficult to pre- 
dict; these are the focus of the present work. 

Recently, various authors have looked into the 
connections between many-body qua ntum dynam- 
ics and classical chaotic behaviorPS In particular, 
it was conjectured in Ref. 8 that microscopic chaos 
plays a role in the long-time behavior of nuclear 
transverse magnetization decays, and that, as a re- 
sult, such decays all exhibit a generic behavior after 
a few times T2, regardless of the initial state of the 
spins at the start of the decay. This long-time be- 



havior has the formPESl 

G{t) ~ e- 7t cos(c^ + 0), (1) 

where G(t) denotes the nuclear magnetization trans- 
verse to the applied magnetic field, <fr is a phase fac- 
tor, 7 is the exponential decay coefficient and u> is 
the beat frequency. Experimentally, this means that 
if the initial transverse spin configuration is varied 
through the application of one or more radio fre- 
quency (rf) pulses with varying phases and inter- 
pulse delay times, then the corresponding transverse 
decays (observed after the final pulse) will have dif- 
ferent initial behavior, but then will eventually all 
converge to the form of Eq. ([!]) with the same values 
of 7 and u>. The parameters 7 and uj have the fol- 
lowing two important properties: 1) each is on the 
order of T^ 1 and 2) each is universal in the sense 
that its respective value is independent of the initial 
transverse spin configuration generated by the rf ex- 
citation pulses. These two properties, together with 
the fact that T2 is the shortest intrinsic timescale 
of the system, imply that Eq. (JlJ cannot be de- 
rived from conventional statistical-mechanical for- 
mulations, which require a separation of timescales 
that does not exist in these systems. The possible 
role of microscopic chaos in the present context is 
discussed in Section Ivl 

Morgan, et alP^ were the first to look for these 
predicted effects with NMR in two solid polycrys- 
talline xenon samples: one composed of 86% 129 Xe 
(isotopically enriched) and the other composed of 
27% 129 Xe (naturally abundant). They found that 
the transverse magnetization decays in these sam- 
ples did indeed approach the universal form given 
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TABLE I: Description of the xenon lattices under 
consideration. They are all face-centered-cubic but 
differ in the concentration of 129 Xe. NMR is 
performed exclusively on the 129 Xe spins in the 
lattice. The xenon isotopes with even numbers of 
nucleons (shown together in the final column) all 
have nuclear spin 1 = 0. 

NanTi System %^ u Xe % ial Xe % even Xe 



43A 


I 


85.6 


1.9 


12.5 


145C 


II 


62.7 


27.3 


10.0 


146A 


III 


54.7 


36.1 


9.2 


146B/145A 


IV 


46.5 


45.2 


8.3 


145B 


V 


29.6 


64.0 


6.4 


123A 


VI 


27.5 


21.5 


51.0 



by Eq. (JlJ by measuring the FID and various solid 
echoes. 

The present work expands the work of Morgan et 
alP3 in the following important ways: 

1) by exploring six different xenon lattice systems 
obtained by systematically varying the ratios of 
xenon isotopes in the samples (see Table [I]), 

2) by exploring the behavior of 19 F NMR in a differ- 
ent material (CaF2) for several orientations of mag- 
netic field, 

3) by employing a wide variety of pulse sequences in 
addition to the solid echo sequences used in Rcf. 11 
to generate distinct initial spin configurations in the 
sample, and 

4) by significantly improving the quality of signals 
over those obtained in previous work. 



II. 



THEORY 



An ensemble of thermally isolated interacting nu- 
clear spins is described by the truncated magnetic 
dipole Hamiltonian in the rotating framed 



G(t) =Tr j^S^(t) 



(4) 



where 



P(t) 



e mt poe~ mt 



(5) 



and po is the density matrix at time t — 0. Here we 
have chosen time origin (t = 0) to be at the time 
corresponding to the end of the final pulse in each 
applied rf sequence. Equation ([I]) is the long-time 
form of Eq. Q predicted by Ref. 8. 

To generate distinct initial density matrices for 
our spin ensembles, we used a variety of pulse 
sequences with varying interpulse delay times r. 
We h ave a nalyzed the FID (90°), various solid 
echoePEl (90° - 

(JB) echo,™j(90° - n 
magic echcPHfil (90 y 



90° 



the Jeener-Broekaert 
45° -T 2 -45° ), and the 



90y - T 



HI 



90„ 



Here the subscripts for the angles indicate the axis of 
the rf pulse- induced rotation in the Larmor reference 
frame, y indicates the negative direction along the 
y-axis. The meaning of H\ is explained at the end of 
Section [H] The solid echo, JB echo, and magic echo 
manipulate spin coherences generated by spin-spin 
interactions and, in this sense, are not like Hahn 
echoes that simply compensate for an inhomoge- 
neous external field. Below we will briefly describe 
each sequence; they are reviewed and analyzed in 
great detail elsewhere.^ 

In the solid echo sequence (90° - t-90°), two 90° 
pulses are applied out of phase with a n interpulse 
delay time r that can be variedP^ESl The density 
matrix which follows the second pulse is given by: 



(6) 



^ \J~krSJ~k,In + ^k^n) + ^kn^-k-^n] > (2) where the "rotated Hamiltonian" Hr is given byf 



.20 



k<n 



where I" is the quantum mechanical operator of the 
n th nuclear spin. The coefficients for the interaction 
of the k th and n th spins are 

T x __ _1 „ „ (l-3cos 2 fl) 
■Jkn- 2 Jkn- 2 |r fe -r„|3 ' W 

where 9 is the angle between the z-axis defined by 
the external static magnetic field and the internu- 
clear vector r-fc — r„. We describe the transverse 
NMR decays by the function G(t) which can be ex- 
pressed as: 



KR = E [ J LiP X n + JUPl + JkJkl'n] (7) 
k<n 

For pairs of interacting spins 1/2, the solid echo 
yields a perfect recovery of the initial magnetization 
at time r after the second pulse, similar to the con- 
ventional Hahn echo. However, in many-spin sys- 
tems complete refocusing does not occur, and the 
deviations from complete refocusing become more 
pronounced as r is lengthened.^ As follows from 
Eq. ([6]), the solid echo pulse sequence creates new 
initial spin density matrices as r is variedP^ 



3 



The JB echo (90° - n - 45° - r 2 - 45° ) par tially 
transforms Zeeman order into dipolar orderl^E^ 
The 45° pulse performs the partial transformation 
into dipolar order — a state similar to that produced 
by adiabatic demagnetization in the rotating frame. 
The 45? pulse generates a JB "echo" from the dipo- 
lar order, better regarded as an FID different from 
the usual one generated from Zeeman order. Here 
we do not explicitly include the expression for the 
density matrix after the final pulse. 

The magic echo sequence is 90 y — r — 90 s — t' h1 — 
90 y . Hi is the amplitude of the radio fre- 



' HI 



quency pulse applied at resonance, and t' h repre- 
sents the duration of time that the pulse is applied. 
Hi represents the amplitude of the resonant pulse 
applied with a 180° phase shift from Hi. We used 
the magic echo as a check on our control of spin 
manipulation; the sequence is designed to reverse 
the Hamiltonian in Eq. |2]) and reproduce the entire 
FID shape. The magic echo is an exceptional pulse 
sequence that completely reverses the spin dynamics 
in the limit of very strong Hi , and therefore the on- 
set of the universal asymptotic behavior occurs not 
after a time on the order of several T 2 , but rather 
after a time t' = t' h — t (peak of the echo) plus 
several times 



III. 



EXPERIMENTAL PROCEDURES 



We explored a wide parameter space in search 
of experimental evidence for the predicted univer- 
sal long-time behavior in isolated quantum systems. 
First, we looked at six different isotopic composi- 
tions of solid hyperpolarized xenon using the solid 
echo sequence with varying values of r. Then, us- 
ing one of these compositions (System I) , we applied 
three different pulse sequences (FID, solid echo, JB 
echo) to observe the long-time behavior associated 
with each initial transverse spin configuration. Next, 
we examined 19 F in CaF 2 single-crystals at three dif- 
ferent orientations of applied magnetic field as well 
as in polycrystalline CaF 2 powder. We used three 
different pulse sequences to observe their long-time 
behavior (FID, solid echo, magic echo). 



A. Solid Xenon 

Hyperpolarized 129 Xe is an ideal system for in- 
vestigating the hypothesis of a universal long-time 
behavior. In this van der Waals solid at 77 K, the 
xenon spin system is thermally well-isolated from its 
environment {T\ T 2 ). The NMR signals from a 
macroscopic sample are large (greater than lmV be- 
fore amplification); the long-time reqime can thus 



be examined with good signal-to-noise-ratio (SNR). 
Additionally, T 2 is long (« 0.5 ms), allowing for rel- 
atively simple radio-frequency manipulation of the 
sample, as the corresponding bandwidths are small 
and the deadtimes short in comparison. These con- 
ditions make the use of short "hard" rf pulses practi- 
cal. Finally the 129 Xe nuclei have spin =1/2, which 
makes them maximally non-classical in terms of spin 
angular momentum quantization, and are therefore 
ideal for exploring quantum connections to conven- 
tionally classical phenomena such as chaos. The 
magnitude of the spin quantum number also means 
we did not have to account for a nuclear quadrupolc 
moment. 

We used dynamic nuclear polarization in the form 
of spin-exchange optical pumping^ to generate nu- 
clear polarizations of wl0% in solid 129 Xe nuclei. 
We then measured the transverse relaxation of these 
nuclei over many orders of magnitude (see Fig. [2| . 
These experiments generate extremely large signals 
that can be precisely measured over unconvention- 
ally long times, i.e., many times T 2 . Spin-exchange 
convection cells^ allowed for rapid and repetitive 
generation of a few millimolcs of highly polarized 
liquid xenon, which we then froze rapidly by expos- 
ing the sample to a stream of liquid nitrogen. During 
rapid freezing, the majority of the polarization in the 
liquid xenon sample survives the phase transition to 
the solidP 

The solid xenon samples were generated in an ap- 
plied magnetic field of 2 T ( 129 Xe Larmor frequency 
24.56 MHz) and maintained at a temperature of 77 
K. Induction signals were acquired with an Apollo 
(Tecmag) NMR spectrometer using 10 //s square ex- 
citation pulses with a receiver dead time of 60 us. 
At 77 K, 129 Xe T x m 2.5 hrs 211 while T 2 0.5 mspl, 
justifying our assumption of thermal isolation dur- 
ing the transverse decay. Our samples consisted of 1 
to 2 mmol xenon in various isotopic ratios (see Ta- 
ble [I]) . We acquired the transverse decay either after 
one 90° pulse (in the case of the FID) or after the 
last of a series of pulses that create various echoes. 
The enormous signals thus acquired were too large to 
fit into the 16-bit digitizer native to the Apollo sys- 
tem. Acquisition was therefore carried out in a two 
separate experiments with two different NMR-gain 
settings. There was sufficient overlap of the tem- 
poral signal in both experiments to ensure that the 
reconstruction of the entire decay was unambiguous. 

To investigate the atomic arrangement of the 
xenon in our samples (see discussion in Section IV 
below), we undertook X-ray diffraction studies on 
a natural xenon sample with the isotopic makeup 
of System VI. A room-temperature sample was ex- 
posed to nitrogen gas at 100 K (well below the 160 
K freezing point of xenon), freezing it suddenly in 
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a manner similar to our NMR experimental proce- 
dure. X-ray diffraction was performed with a 4.1 
kW Bruker Copper K- Alpha X-ray source operating 
at a wavelength of A = 1.5814 A. The sample was 
maintained at 100 K by a built-in cryostat. Typi- 
cal results from one of ten runs are shown in Fig. [I] 
The sharp well-defined rings are characteristic of a 
polycrystalline solid; i.e., randomly oriented crys- 
tals of definite atomic spacing. Analysis of the rings 
due to xenon (after background subtraction) yields 
the known nearest-neighbor spacing of -4.4 AtSlfor 
xenon in a face-centered-cubic (fee) lattice. 

.A 

FIG. 1: (Color online) Result of X-ray diffraction 
experiment of rapidly frozen xenon. Sharp rings 
indicate a polycrystalline sample with well-defined 

lattice spacing, determined to be 4.4 A 
(nearest-neighbor spacing) in agreement with the 
known value P2 



B. Calcium Fluoride 

The CaF 2 crystals used in our experiments were 
obtained from Optovac, Inc (now owned by Corn- 
ing) and are doped with 0.1% yttrium or gadolin- 
ium to reduce T± to « 1 s. These are nominally 
cylindrical samples with the cylinder axis perpen- 
dicular to the [100], [110], or [111] crystal direction; 
these respective directions would thus be parallel to 
the applied magnetic field in a typical transverse 
solcnoidal NMR coil. Because these samples were 
altered somewhat to fit into our NMR probe, the 
crystal directions are no longer exact (though care 
was taken to maintain them as much as possible). 



T 2 25 (j,s for all CaF 2 samples. The signals are 
acquired at 2 T ( 19 F Larmor frequency 83.55 MHz) 
with the same Apollo (Tecmag) NMR spectrometer 
using 2 /is square pulses with a receiver dead time of 
15 (is at room temperature. The CaF 2 powder was 
obtained from J. T. Baker Chemical Co. and con- 
tains nominally 98% CaF 2 with traces of chlorides, 
sulfates, heavy metals, and iron. T\ of the powder 
was measured by saturation-recovery to be ~ Is. 

IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 

A. Solid Hyperpolarized 129 Xe 

1. The Free Induction Decay 




Time (ms) 

FIG. 2: (Color online) 129 Xe FID corresponding to 

System I (86% 129 Xe; see Table [!} shown on a 
semi-log plot. The solid line is a fit of the absolute 
value of Eq. ([I]) to the long-time signal at 
t > 2.5 ms. The inset shows the FID absorption 

signal on a linear scale, illustrating the 
zero-crossings of the FID. Later portions of the 
linear signal have been enhanced as indicated for 
clarity. 

In Fig. [2] we show a typical 129 Xe FID at 2 T and 
77 K. The magnitude of the quadrature detected sig- 
nal is plotted vs. time on a semi-log plot to illustrate 
the decay over many orders of magnitude. The SNR 
here is about an order of magnitude better than the 
previous studyP^ Since the plot shows the log of the 
absolute value of the FID signal, the cusps indicate 
the zero-crossings of the FID (see the inset graph of 
Fig. [2] for a linear-scale plot of the absorption signal 
explicitly showing the zero-crossings). The decay 
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coefficient 7 = 1.25± 0.04 ms" 1 and the oscillation 
frequency uj = 2.06 ± 0.04 rad/ms were extracted 
from a fit of the long-time portion of the FID to 
Eq. ([IJ. These parameters are on the order of T 2 _1 
and agree well with the values obtained by Morgan, 
et. alP^ T 2 was estimated as \J f / Af 2 where M 2 is 
the second moment of the lineshape given by^ 



M, = - 



1 d 2 G{t) 
G(0) dt 2 



(8) 



y=0 



2. The Solid Echo 



In Fig.[3]we show signals associated with ten solid 
echoes obtained for different values of the interpulse 
delay time r together with the free induction decay. 
These data were collected for 86% 129 Xe (System 
I). The initial portion of the FID following the first 
pulse was acquired along with the solid echo follow- 
ing the second pulse in order to provide a relative 
amplitude normalization and thus to correct for po- 
larization levels that varied somewhat from exper- 
iment to experiment. Care was taken to eliminate 
residual longitudinal magnetization from the echo 
signal, and this was checked by using a 90° — r — 90° 
sequence prior to each echo experiment and verify- 
ing a null signal (i.e., precise 90° pulses). We note 
that some of the data reported in Morgan, et alP^ 
from a similar experiment likely are corrupted by a 
small component of magnetization left along the z- 
axis by the first pulse that is then excited into the 
transverse plane by the second pulse. The observed 
7r phase shift of the echoes with respect to the FID 
shown in that work was therefore likely coincidental. 
This does not affect their conclusions concerning the 
universal long-time behavior of the signals, but it is 
critical for exploring relations concerning the ampli- 
tudes and phases of solid echo signals relative to the 
FID in the long-time regimeP^l 

In Fig. [3cJ we make explicit time shifts to the data 
shown in Fig. |3a| and Fig. |3b| to compare the behavior 
of each signal in the long-time regime. The dashed 
line in Fig. 3c is a plot of Eq. ( 1 ) with the parameters 



obtained from the fit in Fig. 2 These values are in 
good agreement with those obtained by Morgan, et 
alP on the same system. Similar long-time fits to 
the individual solid echoes produced values for 7 and 
uj consistent with those obtained for the FID; it is in 
this sense that the long-time behavior is universal, 
in contrast to the varied decay shapes exhibited in 
the initial portions of the decays. 

Tabic |TT] lists the values of the decay constant 7 
and the beat frequency uj for six different xenon iso- 
topic compositions. For each composition, the pa- 
rameters were obtained from a fit to the FID. As 



a check, fits were also performed to each individual 
echo to verify that the extracted parameter values 
were in agreement with those obtained from the fit 
to the FID. 

In Fig. [4] we show the data for System V, repre- 
sentative in SNR of the systems with lower concen- 
trations of 129 Xe. Though the lower isotopic con- 
centration of 129 Xe has reduced the SNR relative to 
System I, we still observe convergence of different 
transverse decays to the respective universal values 
of 7 and uj for each system. The echo amplitudes are 
normalized in the same way as for the other xenon 
compositions as described above. The data similar 
to Fig.[3]and Fig.|4]which pertain to the other xenon 
isotopic concentrations are included in the supple- 
mental section to this workl^l 



TABLE II: Values of the decay coefficient 7 and 
beat frequency uj for each xenon isotopic 
composition tested in our experiments. The fourth 
column refers to the percentage of even (spin-zero) 
xenon isotopes. Data for the compositions marked 
with * are found in the supplemental section [Ref. 
24]. 



System %" y Xe % iJ1 Xe % 


cvon Xe 


7 (ms ') 


uj (rad/ms) 


I 


85.6 


1.9 


12.5 


1.25 ± 0.04 


2.06 ± 0.04 


II 


62.7 


27.3 


10.0 


1.34 ± 0.34 


1.95 ± 0.14 * 


III 


54.7 


36.1 


9.2 


1.42 ± 0.19 


1.96 ± 0.05 * 


IV 


46.5 


45.2 


8.3 


1.24 ± 0.05 


1.71 ± 0.10 * 


V 


29.6 


64.0 


6.4 


1.50 ± 0.06 


1.74 ± 0.05 


VI 


27.5 


21.5 


51.0 


1.03 ± 0.10 


1.53 ± 0.05 * 



3. The Jeener-Broekaert Echo 

Figure[5]shows the JB echo (90°-t 1 -45°-t 2 -45°, 
) for SysTem I along with the corresponding FID 
from Fig. [2] The interpulse delays are t\ = 0.5 ms 
and r 2 — 3 ms. Usually, Ti is much longer than 
T 2 in the JB echo sequence. In our experiments, 
r 2 ~ 6 x T 2 — sufficiently long to generate a spin 
state distinct from that generated by the FID or 
solid echoes. While the initial portion of the JB 
echo is distinct from that of both the FID and the 
various solid echoes, its long-time behavior exhibits 
the form of Eq. ([I]) and with the same values of 7 
and uj (see Fig. [5]) . 
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FIG. 3: (Color online) 129 Xe solid echoes in xenon enriched to 86% 129 Xe (System I). In (a) and (b) we 
show ten solid echoes on a semilog plot together with the FID of Fig. [2j Signals are split between (a) and 

(b) for visual clarity only. In (c), we show the data in both (a) and (b) together with the echoes 
time-shifted to illustrate the convergence of the long-time behavior. The broken line is the fit to the FID, 
also shown in Fig. [2j demonstrating that the solid echoes all exhibit the same long-time behavior as the 

FID with the same values of 7 and u> (see Table |n| . 



B. Calcium Fluoride 

1. The Free Induction Decay and Solid Echoes 

In single-crystal CaF2, different orientations of the 
applied magnetic field H with respect to the crys- 



tal planes imply different interaction coefficients in 
the interaction Hamiltonian T-L [Eq. pi]. Hence, 
experiments performed with each orientation of Hq 
amount to independent checks of the validity of the 
long-time equation Eq. ([!]). In Fig. [6] we show 
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Time (ms) Time(ms) 

(a) (b) 



FIG. 4: (Color online) 129 Xe free-induction decay and solid echoes in xenon having 29.6% 129 Xe (System 
V). Due to lower concentrations of 129 Xe in this system, the SNR is reduced compared to System I. (a) 
Solid echo signals shown on a semilog plot with the FID. (b) The same data as in (a) are shown again with 
the echoes time-shifted to illustrate the convergence of the long-time behavior. The broken line is a fit to 
the long-time portion of the FID from which 7 and ui have been determined (see Table fll|) . 



10V 



10° 



10* 



10 J 



10 V. 



tor 




4 6 
Time (ms) 



10 



FIG. 5: (Color online) 129 Xe FID (solid black line) 
and Jeener-Broekaert echo (blue points) in solid 
xenon enriched to 85.6% 129 Xe (System I). The 
echo is scaled and time-shifted to show the 
uniformity of the long-time behavior. The 
interpulse delays are T\ = 0.5 ms and r 2 = 3 ms. 



a typical 19 F FID at 2 T and 290 K. As with 
xenon, a semilog plot is used and the cusps indi- 
cate zero-crossings. Solid echo data for single-crystal 
[111] CaF 2 (System IX; see Table III) are shown in 
Fig. |S8| These data are representative in behavior 
and SNR of Systems VII - IX. The beat frequen- 
cies in each crystal orientation were first measured 
by Lowe and Norberg^ and the n m ore accurately 
by Engelsberg and LoweP Figure S8 shows that the 
long-time behavior is universal according to Eq. |TJ) 
for 19 F in a single-crystal sample of CaF 2 . The same 
conclusion was drawn from the data for the other two 
crystal orientations! 2 ^ The beat frequency uj mea- 
sured by us for each orientation of Hq is compared 
to that obtained by Engelsberg and Lowe^ in Ta- 
ble |III| The small discrepancies between the values 
obtained by us and those obtained by Engelsberg 
and Lowe are likely due to slight misalignments of 
the samples in our experiment; see Figs. S5c, S6c, 
and S7d in the supplemental section (Ref. 24) for a 
graphical comparison of our data with those of En- 
gelsberg and Lowe. 

Fig. [8] shows the FID and four solid echo signals 
associated with the CaF 2 powder (System X; see 
Table [ill]) . The CaF 2 powder samples (as with all of 
the xenon samples) have many crystallites with ran- 
dom orientations; the transverse decays are there- 
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TABLE III: Values of the decay coefficient 7 and beat frequency 10 for each CaF 2 system tested. Columns 
4 and 5 are our experimental data for the decays. In the last column we list the values of uj calculated from 
the data obtained by Engelsberg and Lowe (taken from Ref. 7) for purposes of comparison. The impurities 
for the powder sample (System X) are as follows: 0.9% CI, 0.5% SO4, 0.3% heavy metals (such as Pb), 
0.3% Fe. Data for the crystals marked with * are found in the supplemental section [Ref. 24]. 



System Direction Impurtiy 7 (ms ) uj (rad/ms) uiel (rad/ms) 

~^VTI [TOO] 0.01 % Gd 51.0 ± 0.1 144.4 ± 0.1 153.32 ± 0.09 : 

VIII [110] 0.01 % Y 43.3 ± 0.1 106.2 ± 0.1 101.9 ± 0.7 * 

IX [111] 0.01 % Gd 31.9 ± 0.1 66.2 ± 0.1 65.6 ± 0.3 * 

X powder 2% n/a n/a n/a 
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FIG. 6: Normalized 19 F FID in CaF 2 
corresponding to the applied field Ho along the 
[110] direction (System VIII; see Table [Til]). The 
solid line is a fit of the absolute value of Eq. (JlJ to 
the long-time signal at t > 45 fj,s. T 2 is approx 25 
[is in this system. 



fore independent of sample orientation. However, 
unlike in xenon, the initial beating seen in the sig- 
nal washes out such that a single decay parameter 
and beat frequency are not present; therefore the 
long-time regime cannot be characterized by a uni- 
versal behavior within our experimental resolution. 
Furthermore, the time-shifts and amplitude manip- 
ulation of the solid echo signals shown in Fig. 8b 
do not give conclusive evidence of a convergence of 
these signals to the behavior of the FID after sev- 
eral times T 2 . Polycrystalline CaF 2 was the only 
system we examined that, within our instrumental 
resolution, did not exhibit the universal long-time 
behavior predicted in Ref. 8. 

Solid xenon samples frozen from the liquid phase 
should also be in a polycrystalline state with ran- 
domly oriented crystalline regions, similar to a poly- 
crystalline powder. Previous work has showrP (and 
we have verified here) that distinct orientations of 



single-crystal CaF 2 with respect to the external 
magnetic field produce different beat frequencies in 
the transverse decay. The frequency differences are 
large enough that mixing them (as is the case in a 
polycrystalline powder) may cause the Lowe beating 
to wash out after a few times T 2 (see Fig. [8| . One 
might have anticipated a similar wash out effect in 
polycrystalline xenon. Instead, the solid xenon sig- 
nals exhibit a very precise beat pattern that lasts 
over several orders of magnitude (see Fig. [2|. This 
phenomenon could be explained if the beat frequen- 
cies of a xenon single-crystal in different orientations 
with respect to the external magnetic field were simi- 
lar or identical, or if the xenon were frozen in a glassy 
state with less order than a rigid polycrystalline 
structure. Indeed, despite some claims in the early 
literature that it should be possible,^ all recent at- 
tempts to vitrify monatomic liquids (noble gases, 
metals, or otherwise) have been unsuccessful to date 
despite achieving freezing rates on the order of 10 10 
K/s and pressures on the order of 10 GPa. The only 
reported exception is a study which achieved vit- 
rification of monatomic germanium by exploiting a 
particular property of tetrahedral liquids that does 
not exist in noble gasesP^ No other monatomic liq- 
uid has yet been vitrified. Our own X-ray diffraction 
data confirm that the frozen xenon is polycrystalline 
with no glassy component (see Fig. [T]) . 

Why the xenon samples exhibit such well-defined 
beats and universal long-time behavior is an out- 
standing question; it is likely related to the differ- 
ence between the fee lattice of solid xenon and the 
simple cubic lattice formed by 19 F in CaF 2 . Each 
lattice site has 12 nearest neighbors in the fee lat- 
tice and 6 nearest neighbors in the simple cubic lat- 
tice, i.e., the arrangement of nearest neighbors is 
more spherically symmetric for the fee lattice. As 
a result, the powder average of the long-time beats 
should be less washed out in the fee case. On the 
other hand, the CaF 2 powder may still exhibit well- 
defined beats at later times (not accessible in this 
experiment) once the crystal grains corresponding 
to the smallest value of gamma begin to dominate 
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FIG. 7: (Color online) 19 F FID and solid echoes in single-crystal [111] CaF 2 (System IX; see Table III I. In 



(a) and (b) we show ten solid echoes on a semilog plot together with the FID. Signals are split between (a) 
and (b) for visual clarity only. In (c), we show the same data as in (a) and (b) with the echoes time-shifted 

to illustrate the convergence of the long-time behavior. 



the overall signal. We pla n to address this question 
theoretically elsewhere.^ Unfortunately, noble-gas 
single-crystals can be generated only with consider- 
able difficulty and in very small sizes™ Furthermore, 
producing single-crystal hyperpolarized 129 Xe would 
appear close to impossible, making direct experi- 
mental NMR measurements of single-crystal solid 
xenon prohibitively difficult. 



2. The Magic Echo 

To contrast our use of the solid echo and JB echo 
to generate distinct initial transverse spin configura- 
tions that then evolve to a universal long-time be- 
havior, we used the magic echo^ (90 a — t — 90^ — 
t' h1 — t^ 1 — 90 y ) in [111] single-crystal CaF 2 as a form 
of control experiment to show that it is possible in 
these systems to reverse the sign of the Hamiltonian 
and reproduce the entire FID. The resulting signal 
from the magic echo sequence shown in Fig. [9] is in- 
deed found to have nearly the same shape as FID 
generated from a single 90° pulse, notably even in 
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FIG. 8: (Color online) F signals in CaF2 powder (System X; see Table III). In (a), solid echoes for four 



different values of r are shown with the FID (black). In (b), the solid echo signals are time- and 
amplitude-shifted to attempt to match the long-time portions of the signals. The Lowe beats wash out 
after several times T 2 , and they are neither similar to each other nor can they be characterized by a single 
beat frequency and decay coefficient within experimental resolution. 



the early portion of the FID decay. For our exper- 
iment, t = 11 /is and t' Hi = r'g = 180 /is. The 
left-hand side of the magic echo near the time ori- 
gin is not perfectly symmetric with the signal at 
t > 200/iS due to residual effects of the strong rf 
excitation pulse. As expected,^ the full dipolar re- 
focusing (and therefore the maximum of the magic 
echo signal) occurs at a time t = t' h — r. 



V. ROLE OF MICROSCOPIC CHAOS 

It has been argued in Refs. [29, 30] that the oc- 
currence of the well-defined long-time behavior (1) 
in the absence of the separation of time scales is 
already a strong indication of the role of the micro- 
scopic chaos. However, such a behavior could also 
appear in approximate theories of spin dynamics as a 
consequence of rather crude approximations. It was 
the prediction of the identical frequency and decay 
constants for the FID and all sorts of echoes^ that set 
apart the theory based on the notion of mi croscop ic 
chaos from other theoretical descriptions l 6 l 10 l 3 IE] 
which, to the best of our knowledge, are not in posi- 
tion to make such a prediction. Below we elaborate 
on these issues. 

Although the long-time behavior (1) is reminis- 
cent of damped harmonic oscillator, one has to ap- 
preciate that damping in the latter case would orig- 
inate from friction, which is usually caused by mi- 
croscopic motion, which, in turn, is much faster than 
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FIG. 9: (Color online) 19 F FID and magic echo in 

single-crystal [111] CaF2 (System IX; see 
Table |!n|. Black line shows the FID, while the blue 
line maps out the signal obtained following the 

magic echo pulse sequence. The time origin 
corresponds to the end of the final pulse in the 
magic echo sequence. The magic echo signal begins 
after a 50/is receiver recovery time. The FID has 
been time-shifted to lie over the portion of the 
magic echo that reproduces the FID. Note that the 
decays are indistinguishable even in the early 
(non-universal) portion of the FID signal. 
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the period and the decay time of the damped oscil- 
lator. In the context of NMR, similar well-defined 
exponential damping appears in the proble ms o f ex- 
change narrowing and motional narrowingPE^l 

It should be further noted that the occurrence of 
well-defined frequency beats in the tail of FIDs and 
echoes is a somewhat counterintuitive effect, which 
is determined primarily by the negative relative sign 
of the J z and J- 1 in Eq. (3). These beats, in par- 
ticular, have little connectio n to the oscillations as- 
sociated with a Pake doubletP^H — the latter orig- 
inates from a highly discrete quantum structure of 
the FID problem for two isolated spins 1/2. In solid 
xenon and in CaF 2 , each spin 1/2 is coupled to the 
whole lattice through many interacting neighbors. If 
the coupling constants J z and J 1 - were of the same 
sign, the FIDs in CaF 2 and solid xenon would likely 
exhibit no beats, while the beats associated with 
a Pake doublet would remain with a modified fre- 
quency value. The above emphasis on correlation 
effects is corroborated by simulations for a lattice of 
classical spins^l where similar beats have been ob- 
served — clearly without any role for quantum dis- 
creteness. These and other factors determining the 
occurrence of beats are discussed in detail in Ref. 30. 
Another difference between the Pake-doublet signal 
and the FIDs in CaF 2 or solid xenon is that the cor- 
responding power spectrum in the former case has 
a well pronounced peak at the oscillation frequency, 
while in the latter case it exhibits no detectable sig- 
nature of the long-time beatsP^ 

In general, FID tails oscillate and decay on the 
fastest natural time scale of the nuclear spin system, 
which is the timescale of individual spin motion in 
the local field created by the neighbors. Historically, 
the absence of the separation of time scales was the 
reason why no universal shape for the FID tails was 
expected a ■priori. Noteworthy in this respect was 
the proposal of Abragam to approximate FIDs in 

_ a 2 t 2 

CaF 2 by function e = sm(bt)/bt, where the pa- 
rameters a and b were to be determined from fitting 
of the exact values of the second and the fourth mo- 
ments P 

The theoretical problem of calculating the FID in 
solids is non-perturbative. Many attempts to deal 
with it resulted in various uncontrollable approxima- 
tion schemes. These schemes aimed at reproducing 
the entire behavior of FID, of which the long-time 
tail was a rather insignificant part. At the same 
time, relatively crude approximations^!^!] made to 
predict the initial FID behavior sometimes resulted 
in the long-time behavior having form (1). 

A particularly strong claim of predicting the long- 
time behavior (1) was made by Borkmans and Wal- 
graef. 32 This claim, however, was criticized as incon- 
clusive in Chapter 8 of Ref. 40, and in the appendix 



of Ref. 30. Borkmans and Walgraef described the 
long-time behavior of FIDs by a memory-function- 
like integro-differential equation with memory func- 
tion (integration kernel) having known initial behav- 
ior but unknown long-time behavior. Looking at the 
initial behavior, they assumed that the entire mem- 
ory function had a Gaussian shape, and thereby ob- 
tained the long-time behavior of form (1). However, 
in fact, the long-time behavior of the FID was en- 
tirely controlled by the unknown long-time behavior 
of the memory function: for the long-time FID be- 
havior (1) to follow from the memory function equa- 
tion, the asymptotic decay of the memory function 
should be no slower than exponential. If it were 
slower than exponential, e.g. power law, then the 
asymptotic FID behavior would also exhibit a power 
law decay rather than (1). 

To the best of our knowledge, the theory of Bork- 
mans and Walgraef has never been applied to solid 
echoes and other pulse sequences investigated in this 
work, and, in fact, we are not sure whether, if ap- 
plied, it would be manageable at all. If it were man- 
ageable, it is further difficult to see how that the- 
ory would reproduce the main result of this work, 
namely, the same long-time constants for different 
pulse sequences. Different pulse sequences would re- 
sult in different initial conditions and hence differ- 
ent parameters for the memory functions. It then 
appears rather unlikely that different memory func- 
tions would lead to the same long-time decay con- 
stants. A similar critique is applicable to all other 
approximations attempting to predict the long-time 
behavior of NMR signals from the knowledge of a 
few initial time derivatives. 

The prediction of quantitatively identical long- 
time behavior for different pulse sequences was not 
made until the long-time behavior of NMR signals 
was conside red o n the basis of the notion of micro- 
scopic chaosP^ The analysis of Ref. 30 has focused 
exclusively on the fundamental reasons behind the 
fact that the long-time functional form (1) is not 
a crude approximation but rather a robust accu- 
rate property. This analysis predicted that non- 
integrable lattices of interacting classical or quan- 
tum spins universally exhibit the long-time behavior 
(1), while the integrable cases amenable to the an- 
alytical calculations of FIDs exhibit non-universal 
behavior different from (1). Non-intcgrability in 
the interacting spin systems is the rule, while in- 
tegrability is an exception. In classical spin sys- 
tems, non-integrability is generally associated with 
microscopic chaos caused by non-linear interaction 
between spins. Mathematically, chaos is defined as 
an exponential instability of the phase space trajec- 
tories with respect to small differences of the initial 
conditions. Calculation of classical FIDs requires 
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averaging over all possible chaotic phase space tra- 
jectories. In the quantum case, the notion of chaos 
cannot be defined at the level of exponential insta- 
bilities in the phase space because the Heisenberg 
uncertainty principle precludes one from defining a 
phase space trajectory. However, the connection be- 
tween classical and quantum chaos can be conjec- 
tured for spin lattices at the level of spin correlation 
functions.™ 

According to the analysis of Ref. 30, almost any 
initial probability distribution in a chaotic classical 
spin system would evolve to exhibit universal pat- 
terns associated with stable and unstable directions 
in the phase space. These patterns are complicated 
and intrinsic to a given Hamiltonian. They exhibit a 
dynamically developing small parameter represent- 
ing the ratio of characteristic scales along stable and 
unstable phase space directions. This small param- 
eter, in effect, substitutes for the separation of time 
scales and leads to Markovian long-time behavior 
of ensemble-averaged quantities on non-Markovian 
timescales [e.g. Eq.(l)]. As time progresses, the 
above patterns become increasingly singular and are 
expected to underlie chaotic eigenmodes of the time 
evolution operator. These e igenm odes are known 
as Pollicott-Ruelle resonance J 41 * 42 1 Their singular 
character is masked by ensemble averaging, but what 
remains after averaging are the exponential decays 
with or without oscillations. In terms of spin corre- 
lation functions, the formation of the above patterns 
imply the transfer of weight from lower-order to 
higher-order correlations (spin coherences in NMR 
language). It was conjectured in Refs. [8, 20, 30] that 
the same transfer of spin correlations also occurs in 
quantum spins systems and underlies the observed 
exponential decays. The slowest of these decays then 
controls the long-time behavior in a given symmetry 
class of initial conditions. (NMR signals studied in 
this work correspond to the initial conditions hav- 
ing zero wave vector symmetry, i.e., all spins are 
equivalent from the viewpoint of initial conditions.) 
It was estimated in Ref. 30 that the onset of this 
single exponential behavior occurs in the NMR con- 
text on the scale of a few times T-x- The predic- 
tion of the identical long-time behavior for different 
pulse sequences was based on the fact that the above 
Pollicott-Ruelle resonances and the associated expo- 
nential decays are intrinsic for a given Hamiltonian 
but not to the initial conditions, and therefore the 
exponential decay associated with the slowest res- 
onance would dominate the behavior of many-spin 
density matrices and hence would manifest itself in 
the behavior of all observable quantities in the same 
symmetry class independently of the initial condi- 
tionP 

At present, the theory of chaos in both classi- 



cal and quantum many-body systems is still at an 
early development stage; the use of conjectures in 
the above analysis is therefore unavoidable. How- 
ever, since the conjectures are not rigorously proven, 
it is essential to substantiate theoretical claims by 
extensive experimental tests, and this is what the 
present work does. We also hope that the robust- 
ness of the long-time behavior (1) would stimulate 
the search for the theoretical methods of controllable 
calculations of parameters 7 and u> in Eq.(l). These 
methods, if found, are likely to have impact much 
beyond NMR theory. 



VI. CONCLUSION 

The present paper reports on a comprehensive in- 
vestigation of the long-time behavior of NMR sig- 
nals well beyond the regimes explored in previous 
work, thereby establishing a new frontier of what 
is known on this subject experimentally. In com- 
parison with the previous studyp^ which investi- 
gated two polycrystalline samples of one material 
(solid xenon) in a smaller signal intensity range, the 
present study includes polycrystallites/powders and 
single crystals of two different materials (solid xenon 
and CaF 2 ) with a better signal resolution and with 
a much broader parameter space (multiple isotopic 
xenon concentrations, varied external fields for CaF2 
crystals, and a wide variety of pulsed spin manip- 
ulations). In all cases except for the CaF 2 pow- 
der, we have observed sinusoidally modulated expo- 
nential long-time decays characterized by the same 
time constants irrespective of distinct initial condi- 
tions imposed by various spin manupulations. In 
the CaF2 powder, the sinusoidal Lowe beat pattern 
of 1 F washes out within our experimental resolu- 
tion and therefore cannot be characterized by Eq. 
(1). In principle, initial washing out of Lowe beats 
is expected in a powder where the beat frequencies 
of the various crystallites are distinct. In light of 
these complications associated with powder averag- 
ing, the fact that the well-defined universal behavior 
(1) is seen in the polycrystalline xenon systems is 
intriguing and needs to be better understood. But 
more importantly, the results on single crystal CaF 2 
first reported here emerge as more conclusive evi- 
dence of the identical long-time behavior of FIDs 
and various echoes. Since it is also not known a 
priori how far the exponential-oscillatory long-time 
behavior can extend in time, particularly in pow- 
der samples, the improvement in the signal-to-noise 
ratio allowing one to explore longer time ranges is 
another essential advance of this study. 

The bulk of this work provides evidence in sup- 
port of the conjecture that a universal behavior will 
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dominate the long-time transverse relaxation of ma- 
nipulated spin coherences. This conjecture is based 
on the notion of microscopic chaotic mixing induced 
by the interaction between nuclear spins. We note 
that as our system for hyperpolarizing xenon nu- 
clei continues to improve, the resulting highly spin- 
polarized lattice will generate yet another class of 
unique transverse spin configurations that depart 
from the high-temperature limit; these can then be 
examined for further verification of the universal 
long-time behaviorPS 
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System II (62.7% 129 Xe) 
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FIG. SI: (Color online) 129 Xe free-induction decay and solid echo in xenon having 62.7% 129 Xe (System II; 
see Table I of main article), (a) Solid echo signal shown on a semilog plot with the FID. (b) The same data 
as in (a) are shown again with the echo and FID time-shifted to illustrate the convergence of the long-time 

behavior. 
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System III (54.7% 129 Xe) 
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FIG. S2: (Color online) 129 Xe free-induction decay and solid echoes in xenon having 54.7% 129 Xe (System 
III; see Table I of main article), (a) Solid echo signals shown on a semilog plot with the FID. (b) The same 
data as in (a) are shown again with the echoes time-shifted to illustrate the convergence of the long-time 

behavior. 
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System IV (46.5% 129 Xe) 
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FIG. S3: (Color online) 129 Xe free-induction decay and solid echoes in xenon having 54.7% 129 Xe (System 
IV; see Table I of main article), (a) Solid echo signals shown on a semilog plot with the FID. (b) The same 
data as in (a) are shown again with the echoes time-shifted to illustrate the convergence of the long-time 

behavior. 
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System VI (27.5% 129 Xe) 



FID 

0.56 ms 




Time (ms) 

(a) 



Time (ms) 



(b) 



FIG. S4: (Color online) 129 Xe free-induction decay and solid echoes in xenon having 27.5% 129 Xe (System 
VI; see Table I of main article), (a) Solid echo signal shown on a semilog plot with the FID. (b) The same 
data as in (a) are shown again with the echoes time-shifted to illustrate the convergence of the long-time 

behavior. 
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System VII ([100] CaF 2 crystal) 
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FIG. S5: 19 F FID and solid echoes in single-crystal CaF2 with external field along [100] (System VII; see 
Table III of main article), (a) Solid echo signals shown on a semilog plot with the FID. (b) The same data 
as in (a) are shown again with the echoes time-shifted to illustrate the convergence of the long-time 
behavior, (c) FID obtained in this work plotted with that obtained in Ref. 7. Discrepancies in the beat 
frequency can be attributed to slight differences in crystal alignment. 
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System VIII ([110] CaF 2 crystal) 
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FIG. S6: 19 F FID and solid echoes in single-crystal CaF2 with external field along [110] (System VIII; see 
Table III of main article), (a) Solid echo signals shown on a semilog plot with the FID. (b) The same data 
as in (a) are shown again with the echoes time-shifted to illustrate the convergence of the long-time 
behavior, (c) FID obtained in this work plotted with that obtained in Ref. 7. Discrepancies in the beat 
frequency can be attributed to slight differences in crystal alignment. 
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System IX ([111] CaF 2 crystal) 
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FIG. S8: (Color online) 19 F FID and solid echoes in single-crystal CaF2 with external field along [111] 

(reproduced from Fig. 7 in main atricle). In (a) and (b) we show ten solid echoes on a semilog plot 
together with the FID. Signals are split between (a) and (b) for visual clarity only. In (c), we show the 
same data as in (a) and (b) on a semilog plot with the echoes time-shifted to illustrate the convergence of 
the long-time behavior, (d) FID obtained in this work plotted with that obtained in Ref. 7. Discrepancies 
in the beat frequency can be attributed to slight differences in crystal alignment. 



